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ABSTRACT 

Recent simulations of Population HI star formation have suggested that some fraction form in binary 
systems, in addition to having a characteristic mass of tens of solar masses. The deaths of metal-free 
stars result in the initial chemical enrichment of the universe and the production of the first stellar- 
mass black holes. Here we present a cosmological adaptive mesh refinement simulation of an ovcrdcnsc 
region that forms a few 10^ M0 dark matter halos and over 13,000 Population HI stars by redshift 
15. We find that most halos do not form Population HI stars until they reach A/vir ^ 10^ Mq because 
this biased region is quickly enriched from both Population HI and galaxies, which also produce high 
levels of ultraviolet radiation that suppress H2 formation. Nevertheless, Population HI stars continue 
to form, albeit in more massive halos, at a rate of 1O~''M0 yr~^ Mpc~^ at redshift 15. The most 
massive starless halo has a mass of 7 x 10'' Mq, which could host massive black hole formation through 
the direct gaseous collapse scenario. We show that the multiplicity of the Population HI remnants 
grows with halo mass above 10® Mq, culminating in 50 remnants located in 10^ Mq halos on average. 
This has implications that high mass X-ray binaries and intermediate mass black holes that originate 
from metal-free stars may be abundant in high-redshift galaxies. 

Subject headings: cosmology - methods: numerical - hydrodynamics - radiative transfer - star for- 
mation 



1. INTRODUCTION 

The first generation stars (Population HI) form from 
metal-free gas in dark matter halos with M ~ 10^ Mp and 
have a large characteristic mass fe.g. [Abel et"aI1l200l 



Broinin et al.ll2002t [O^Shea fc NormanI 
20091 : IGreif et all I2012t) 



2007tlTurk et atl 

Due to their high mass. 



tic medium (IGM), halo mergers, and cosmological 
accretion. This topic has been extensively stud- 
ied with semi-analyt i c models (iScannapieco et al.l [20031: 
Yoshida et all |2004| iTuinhnsonI l2006t iSalvadori et al.l 



they have short lifetimes (jSchaered l2002[ ) and may go 
supernova, and enrich their surrounding intergalactic 
medium (IGM). For these metal- free stars, Type II super- 
novae (SNe) happen for an initial mass between 10 and 
40 Mq, and much more energetic pair- instability SNe 
(PISNe) might occ ur in stars between 140 and 260 Mq 
IHeger et al.r(j2003D . Once the metallicity passes some 
critical metallicity, ^ 10~^ Zc^ if dust cooling is efficient 



20071 : iKomiva et al. Il2010n. post-proces sing of numer i- 
cal simulations ([Karlsson et al. 2008; Trc nti et al.ll2009l) . 
and direct numerical simulations (Tornatore et al. '2007" 



Ricotti et al. 2 008: Maio et al,, 2 010; Wise et al 20121 
Muratov etlj] 1201^ . For example , iTrenti et all (|200S 



: coolmg 1 

(jOmukai et al.l [200l [Schneider et"an [2006t IC lark et all 
20081) or ~ 10"^ -^ Zq otherwise (jBromm et al, .2001] : 



Smith et al.|[2009l ). the gas can cool rapidly and lower its 



Jeans mass. These metal-enriched (Population II) stars 
have a lower characteristic mass scale and most likely 
have an initial mass function (IMF) that resembles the 
present-day one. However at high redshift, heating from 
the cosmic microwave background (CMB) may limit the 
radiative cooling, and thus increasing the Jeans mass, re- 
sulting in an IMF that also favor s massive star formation 
(ILarsonll2005t ISmith et al.l[2009l) . 

The transition from Pop HI to Pop II star for- 
mation is solely dependent on the metal enrichment 
from the Pop HI SN remnants in the future star 
forming halos. Metal enrichment involves complex 
interactions between SNe blastwaves, the intergalac- 
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suggested that Pop HI stars may still form at the late 
epoch of z = 6 in the under dense regions of the universe 
by post-processing of co smolog i cal sim ulations with blast 
wave models. Murato v et al.l (|201^ also showed that 
Pop HI stars continue to form until z — 6 using direct 
cosmological simulations. 

In addition to the metal enrichment from Pop HI stars, 
heating and ionizing effects from their radiation are cru- 
cial to modeling early structure formation of the Uni- 
verse. ([Haiman et al.l [20001) . Lyman- Werner photons 
from Pop HI stars photodissociate H2 by the Solomon 
process and suppress the formation of Pop HI stars in 
low-mass halos. The higher energy UV radiation from 
Pop HI stars then affects the subsequent structure forma- 
tion through heating and ionizing the surroundin g IGM 
jMachacek et al. 2001; Yoshida et al, 2003; Wise &: Abell 
[20071: lO'Shea fc Norman! 120081 ). Pop HI more mas- 
sive than 260 Mq or less massive than 140 M q may 
direc t collapse to form black holes iHeger et al.l (BHs; 
l2003f ). Accretion onto these massive P op HI BHs is 
a fea sible way to form z > 7 quasars ([Johnson et al.l 
I2OI2D and is also an important source of X-ray radia- 
tion in high redshifts. X-rays from the accretion onto 
Pop HI BHs jA lvarez ct al. 200§) or from Pop HI bina- 
ries|Turk eTaLII2009t iStacv et aL|[20lol:lStacv fc Bromml 
[Ml " may preheat and pre-ionize a large volume of the 
IGM l|Ricotti fc Ostriker|[20p: IMesinger et~aI][20Tl . 
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Figure 1. Snapshots of the refined regions and the most massive halos at redshifts 25 (left), 17.91 (middle) and 15 (right). The images are 
the density-weighted projections of baryon density in cubic volumes 6.6 comoving Mpc on a side (first row), enclosing the refined regions, 
the density-weighted projections of baryon density (second row), temperature (third row) and metallicity (fourth row) in cubic volumes 
10.0 proper kpc on a side, enclosing the most massive halos. 



The impact of Pop III feedback on cosmic evolution 
is dependent on the properties of their host halos, espe- 
cially their masses ([Whalen et al.l 120081 : IMuratov et al.l 
[20Ta) . The sizes of the host halos determine the dis- 
tance the metal from Pop III SNe reach with their blast 
waves, and the escape faction of their UV radiation. This 
makes a detailed study of the Pop III star and remnant 
distribution over a wide range of high redshift galaxies 
necessary. 

It is impossible to observe of Pop III stars during their 
lifetime at high redshift, but they might still possibly be 
detected directly by looking for their PISN explosions be- 



fore their death. This idea has been studied and s hown 
to be p romising for both LSST (iTrenti et aLll2009D and 
JWST (jHummel et aLllMl IWhalen et al.ll20m Un- 
derstanding the population and distribution of Pop III 
stars and remnants in high redshift galaxies is helpful for 
preparing the observation of these events. 

In this paper, we focus on the the formation of Pop 
III stars and the population and multiplicity of Pop III 
stars and remnants in high redshift galaxies. We have 
performed a simulation of a survey volume of over 100 
comoving Mpc'^ that includes a full primordial chemistry 
network, radiative cooling from metal species, both Pop 
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II and Pop III star formation and their radiative, ther- 
mal, mechanical and chemical feedback. The simulation 
runs on such a large volume from cosmological initial 
conditions, so that we can follow the formation and fate 
of Pop III stars in a statistically significant number of 
halos with a wide range of masses from a few million 
solar mass to one billion solar mass. We first describe 
our simulation model in Section 2. Then in Section 3, 
we present our results of the Pop III stars and remnants 
distribution over early galaxies. We discuss the findings 
and possible bias in our simulation in Section 4. 

2. SIMULATION SETUP 

We perform the simulation using the adaptive mesh re- 
finement (AMR) co smological hydrodynam ics code Enzo 
(lO'Shea et al.ll2004[) . Adaptive ray tracing (jWise &: Abell 
|2011|) is used for the radiation transfer of ionizing radia- 
tion, which is coupled to the hydrodynamics and chem- 
istry in Enzo. The chemistry, cooling, and star formation 
and f eedback models use d in this simulation are the same 
as in lWise et all ()2012b[ ). 

We generate the initial conditions fo r the simula- 
tion using Music (iHahn fc Abell 120 111) at z = 99 
and use the cosmological parameters f rom the 7-year 
WMA P ACDM+SZ+LENS best fit (|Komatsu et al.l 
I20T1I) : 17m==0.266, flA = 0.734, 17b=0.0449, h=0.71, 
(78=0.81, and n=0.963. We use a comoving simulation 
box of (40 Mpc)"^ that has a 512'^ root grid resolution 
and three levels of static nested grids. We first run a 
512^ N-body only simulation to z = 6. Then we select 
the Lagrangian volume around two ~ 3 x halos at 

z = 6, and re-initialize the simulation with 3 more nested 
grids to have an effective resolution of 4096^ and an ef- 
fective dark matter mass resolution of 2.9 x IO^Mq inside 
the highest nested grid of 300 comoving Mpc'^. During 
the course of the simulation, we allow a maximum re- 
finement level / = 12, resulting in a maximal resolution 
of 19 comoving kpc. T he refinement crite ria employed 
are also the same as in I Wise et al.l (j2012bf ). but the re- 
finements higher than the static nested grids are only 
allowed in the region that contains only high resolution 
particle, that is ^ 150 comoving Mpc^ at z = 15. At this 
time, the simulation has a large number (■^lOOO) of halos 
with M > 10® M0, where new formation of Pop III stars 
declines rapidly, for statistical analysis. We use results 
at this redshift for our current study. The simulation has 
1.3 billion computational cells and required more than 10 
million CPU hours, and there are three > 10^ Mq halos 
in the refined regions at this time. We will continue this 
simulation to lower redshift for the study of the Pop II 
stars and high redshift galaxies. 

Both Pop II and Pop III stars form in the simulation, 
and we distinguish them by the total metallicity of the 
densest star forming cell. Pop III stars are formed if 
[Z/H] > —4, and Pop II stars are formed otherwise. We 
use the sa me star formation m odels and most of the pa- 
rameters in I Wise et al.l (|2012bD . as well as feedback mod- 
els. The only difference between the two simulations is 
that we use a characteristic mass of 40 Mq for the Pop 

III IMF, which is more in line with the latest results 
of Pop III forma tion simulations (e.g. iTurk et al.l [20091 : 
iGreif et al.ll2012D . ins tead of 100 M,^. Pl ease see Section 
2.2 and Section 2.3 of IWise et all ()2012bD . for the details 
of the star formation and stellar feedback used in the 
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Figure 2. Top panel: Evolution of the cumulative number of 
Population III stars and remnants and the averaged Pop III star 
number formation rate inside the survey volume. Bottom panel: 
Evolution of the total Pop II star mass, and the Pop II and Pop 
III star formation rate by mass. The mass formation rate of Pop 
III stars depends on the choice of M^har- 

simulation, respectively. 

We show the evolution of the high-resolution region in 
Figure 1 at three redshifts, 25, 17.91, and 15. Here we 
show the large scale structure in the inner 6.6 comoving 
Mpc in the top row, and in the remaining rows, we focus 
on the most massive halo at each redshift by showing 
their density-weighted projections of gas density, tem- 
perature, and metallicity. 

3. RESULTS 

We illustrate the evolution of the number of Pop III 
stars and remnants, as well as their formation rate from 
the birth of the first Pop III star to z = 15 in the 
top panel of Figu re 2. Massiv e Pop III stars have very 
short lifetimes (jSchaered 120021 ) and end their life by ei- 
ther direc tly collapsing to bl ack holes or exploding as su- 
pernovae (|Heger et al.ll2003[ ). depending on their initial 
masses. More specifically, they die as Type II SNe if 11 
< M^ / Mq < 40 and as PISNe if 140 < M^ / Mq < 260, 
where M^ is the initial stellar mass, or become black holes 
if their masses are not in these mass ranges. Through- 
out the paper, we use "Pop III remnants" to refer to all 
remains of dead Pop III stars, regardless whether they 
become black holes or supernova remnants, and we use 
"Pop III" to refer to both living stars and dead remnants. 
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Figure 3. Left panel: Number of Pop III, total halos and lialos with Pop III as functions of halo mass at redshifts 17.91 (Top) and 15 
(Bottom). The number of living Pop III stars are stacked over the Pop III remnants showing in the same bin. Right panel: Fraction of 
halos hosting Pop III and average number of Pop III per halo at redshift 17.91 (Top) and 15 (Bottom). 



In the case that remnants have neghgible masses after SN 
events, their star masses are replaced with a very small 
mass proportional to their initial masses and the star 
particles are kept in the simulation to follow their rem- 
nant kinematic distribution. The first Pop III star forms 
at redshift z ^ 29.7. Rates of Pop III star formation 
steadily increase from lO"'^ to above 10~^ stars per 
year per comoving Mpc"^. The Pop III star formation 
rate shows some signs of saturation at z ~ 15 in this 
overdense region of the universe. K\, z — 15, the entire 
survey volume of 138 comoving Mpc'^ contains 13,123 
Pop III stars and remnants and 7,677 halos more mas- 
sive than 5 x 10^ M0. The number densities of Pop III 
and halos that host Pop III are 95 and 55 per comoving 
Mpc'^, respectively. 

The total mass of Pop II stars and the SFR of Pop II 
and Pop III arc shown in the bottom panel of Figure 2. 
There is '--^ 3 x 10® Mq mass in Pop II stars at z = 15. 
The star formation histories f or both Pop II and Pop 
III stars arc similar to those in lWise et al.l ()2012bD . but 
shifted to higher redshifts. We will study the details of 
Pop II star formation in this simulation in a forthcoming 



paper. 

We show the number of Pop III, all halos, and star- 
hosting halos as function of halo mass inside the survey 
volume at z = 17.91 and z = 15 in the left panels in 
Figure 3, while the numbers of living Pop III stars are 
stacked over remnants in the same bin. In the right pan- 
els in Figure 3, shown are the fraction of halos hosting 
Pop III and the average Pop III per halo (over halos hav- 
ing Pop III) as functions of halo mass. 

As our survey volume is chosen over a high density re- 
gion, the halo number density is higher than the Universe 
mean. The halo mass functions fro m our simulat i on are 
well fit with the fitting function of I Warren et al.l ()2006l ) 
boosted by a factor of five for halos heavier than 5 x 
10^ Mq in both redshifts. The number of halos drops 
for halos < 5 x 10^ Mq, due to the lack of dark matter 
mass resolution. At z = 15, our simulation has a factor 
of a few more halos than the prediction from the fitting 
function for halos more massive than 1 x 10® A/q, show- 
ing that the baryon becomes important in massive halos. 
There arc 3 halos with masses over 10^ Mq at z = 15, 
while the first one appears at redshift z ~ 15.8. 
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Figure 4. Phase plots of halo number indicating the distribution 
of Pop III stars and remnants over halos as a function of PopIII 
number and their host halo mass at redshifts 17.91 (Top) and 15 
(Bottom). 

The Pop III distribution over halo mass shows httle 
evolution, except that there are more Pop III appearing 
in high mass halos as those halos building up with time. 
The number of Pop III peaks at ^ 3 x 10^ M©. The 
fraction of halos hosting Pop III also has little change be- 
tween redshifts 17.91 and 15 because the chances to host 
Pop III are very small for halos less massive than 10^ 
Mq in this biased region. Then they gradually increases 
with halo mass and reach 100% for halos ^ 10® M©. In 
this particular mass range of 10^ — lO^M©, halos can cool 
and collapse regardless of the LW background and form 
stars (e.g. IWise & Abef .2007; O'Sh ea & Norman 200S; 
iShang et allloT^Wolcott-Green et al.ll201l[ ). The ex- 
act timing of Pop III star formation depends on the lo- 
cal strength of the LW intensity, causing the gradual in- 
crease in the fraction shown in Figure 3. In addition 
to H2 suppression, nearby galaxies and Pop III super- 
nova can chemical enrich halos so that they form metal- 
enriched stars in their first star formation event without 
hosting any prior Pop III star formation. Afterwards, as 
halos grow through hierarchical mergers, the number of 
Pop III per halos increase with both redshift and halo 
mass. At 2; = 15, halos between 10'' Mq and 10* Mq 
have an average of 10 Pop III, and the number rises to 
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Figure 5. Phase plots of Pop III number indicating the relation 
between Pop III formation and host halo mass as a function of 
PopIII star formation time and their host halo mass at redshifts 
17.91 (Top) and 15 (Bottom). 

about 50 for halos 10^ Mq. 

Shown in Figure 4 are the number distribution of halos 
as a function of halo mass and number of Pop III stars 
and remnants per halo. The number of PopIII that is 
hosted by halos inside a mass bin is very scattered, but 
there is a clear trend that more massive halos host more 
Pop III stars and remnants. The same trend is shown 
in the average number of Pop III per halo in Figure 3. 
At redshift z = 15, the most massive halo without any 
Pop III stars or remnants is 7.16 x 10^ Mq. On the low 
mass end; the least massive halo hosting Pop III is only 
1.82 X 10^ Mq. But this Pop III star forms at about 50 
Myr ago, so this star most likely forms in another more 
massive halo, then is stripped out with its current host 
halo from the massive one. 

In Figure 5, we show the number distribution of Pop 
III as function of their host halo mass and their forma- 
tion time. As old Pop III remnants can be found in halos 
of all masses, young stars and remnants can only be seen 
in low mass halos. No Pop III star forms in halos more 
massive than 3 x 10* Mq. The gas in massive halos 
is enriched by their earlier Pop III, so the further for- 
mation of stars transitions to Pop II. The least massive 
halos having newly formed Pop III stars is ~ 4 x 10^ 
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Mq, as the Pop III formation in smaller halos is sup- 
pressed by the LW background. But our results for Pop 
III formation in low mass halos might be biased. Due to 
the relatively low dark matter (2.9 x 10'^ M©), we may 
miss the formation of M < lO^M© halos in our simula- 
tion. We discuss this issue in the next section and show 
that very little Pop III star formation is missed by the 
current simulation due to resolution effects. 

4. DISCUSSIONS AND CONCLUSIONS 

We present results of a cosmological simulation repre- 
senting a high density region of an unprecedented large 
volume in the early universe to study the Pop III star and 
remnant population and multiplicity within high redshift 
galaxies. The highly resolved simulation volume covers 
thousand of halos and Pop III stars and remnants, al- 
lowing us to make a detailed study of the Pop III star 
formation history and the Pop III stars and remnants 
distribution over a wide range of halo masses. We ob- 
serve a continuous Pop III star formation from z ^ 30 
to z 15, when the calculation was paused for analysis. 
We expect Pop III star formation to continue to lower 
redshifts, but this needs to be verified by running the 
simulation further. We find that the number of the Pop 
III stars and remnants peaks in halos with masses of a 
few X 10^ Mq during this time, while the ratio to host 
Pop III and multiplicity of Pop III increase with the halo 
mass. The Pop III inside a massive halos are likely com- 
ing from mergers of small ones, instead of forming in the 
big ones. New Pop III stars only form in halos between 
4 X 10^ to 3 X 10® Mq, while the formation rate peaks 
at ~ 10^ Mq. So even though there are tens of Pop III 
remnants in a halo of 10^ M©, little to no feedback might 
be expected coming from their still living Pop III stars. 

The dark matter mass resolution in this simulation is 
not high enough to capture all of the Pop III star forming 
halos smaller that a few xlO^ M©, so we could be un- 
derestimating the number of Pop III in our simulation. 
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To inve stigate this, we use the simulation in lWise et al.l 
()2012al) to estimate the missing Pop III stars and rem- 
nants in low mass halos. We show the histogram of halo 
mass in that 1 comoving Mpc^ box simulation at z ~ 7.28 
in Figure 6. Due to the much higher dark matter mass 
resolution, there are many more halos at a few 10^ M©, 
extending to 2 x 10^ Mq. Since the simulation box 
is small, the halo m ass function is not well fit by the 
iWarren et al.l (|2006f ) function by having too many low 
mass halos and too few high mass halos, but still has the 
number of halos ~ 10^ Mq close to cosmic mean. Though 
there are much less halos and Pop III in this smaller simu- 
lated volume, it has a similar distribution of Pop III stars 
and remnants as in our current simulation. Halos with 
-^'jfhaio ^ 2 x IO^Mq do not host any Pop III stars or rem- 
nants in this small box simulation. Pop III star forma- 
tion in these low-mass halos is suppressed when the II2 is 
dissociated by LW radiatio n from nearby stellar sources 
and the background ( Macha cjk et al.ll2001l: iWise &: Abel 
120071: lO'Shea fc Norniai], .20081 ) . The first star forms at 
z ~ 20, when the LW background is already at an inten- 
sity of ~ 10"^'^ erg s~^ cm~^ Hz~^ sr~^, and suppresses 
any H2 formation and thus Pop III star formation in 
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Figure 6. Number of Pop III stars and remnants, all halos, and 
halos with Pop III as functions of hal o mass of the 1 comoving 
Mpc-^ simulation in IWise et"aLl I l2012bl ') at redshift 7.28. With a 
dark matter mass resolution that is finer by a factor of 16, the halo 
population extends from 10^ Mq to 2 x 10^ Mq. However, there 
are no Pop III present in those small halos, which are suppressed by 
a LW background. Indeed, the distribution of Pop III is similar to 
our current results shown in Figure 3, though two survey volumes 
are in very different density fluctuations. 



Figure 7. Top panel: Evolution of the volume averaged Lyman- 
Werner radiation intensity from stars in the simulation. Bottom 
panel: Complementary cumulative distributions of the LW inten- 
sity at 2 = 15. J21 is the LW intensity in the unit of 10~^^ erg s~^ 
cm~^ Hz~^ sr~^. The averaged LW intensity gradually increases 
to and maintains above 10~^ erg s~^ cm"'^ Hz~^ sr~^ for the 
last 40 Myr. This LW radiation is not enough to significantly slow 
down the Pop III formation. 
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halos below a mass of IO^Mq. At later times, we also 
confirm that the minim um star- forming ha lo mass with 
this LW background in I Wise et all (1201 2af ) at any time 
is > 1O^M0. Thus, we conclude that the Pop III for- 
mation rate is small in halos below this mass scale and 
that we do not miss a significant fraction of Pop III star 
formation in our simulation because of mass resolution. 

In this large volume simulation, only the LW flux from 
point sources is included, while the LW background is ig- 
nored. We show the evolution of averaged LW intensity 
and the fraction distribution of LW intensity at z = 15 
over the refined volume in Figure [71 The averaged LW 
radiation intensity builds up gradually and reaches 10~^^ 
erg s""'^ cm~^ Hz~^ sr~^at redshift z ^ 17, and maintains 
at this level for more than 40 Myr. At z = 15, the LW 
intensity in ^ 70% of the volume is higher than 10~^^ erg 
s~^ cm~^ Hz"^ sr~^, while only less than 0.01% of the 
volume (0.014 comoving (Mpc)'^) is filled with LW radi- 



ation stronger than 10 erg s 



-1 sr-^ This 



LW intensity, locally produced even in a rare high density 
region, is still weaker than the expected LW background 
(iWise fc Abell [20051 : iWise et all [20T2al ). and, though it 
suppresses Pop III formation in 10^ Mq halos, is not 
strong enough to have significant impact on the Pop III 
formation in larger halos. 

An important di fi'erence between the simulation in this 
paper and that of IWise et al.l (|2012al) is that the char- 
acteristic mass of Pop III stars here is much lower at 
Mchar = 40 Mq, instead of Mchar = 100 Mq. The new 
choice of characteristic mass results that much fewer Pop 
III stars have initial mass between 140 and 260 Mq, and 
end their life by PISN. This then significantly reduces the 
metal generated by Pop III stars. For exa mple, a 40 Mq 
hype rnova produces only 8.6 Mq metals (iNomoto et al.l 
|2006[) . con iparing to 85 Mq of me tals generated by 180 
Mq PISN (jkeger fc Wooslevl[200l . It is likely due to the 
lower level of metal generation, that Pop III stars con- 
tinue to form in the larger halos (> 10^ Mq) in current 
simulation. On the contrary, there is no Pop III formed 
in halo s more massive than 5 x 10^ Mq in [Wise et all 
(l20T2bl) . 

Because our simulation covers a rare high density re- 
gion, our results are not representative of the cosmic 
mean. We fit the halo mass f unction with the fitting 
function of I Warren et al.l (j2006D . It shows that the halo 
number density at z = 15 is about five times that of 
the cosmic mean at the same redshift or about the cos- 
mic mean at z ~ 10, while the halo number density at 
z = 17.91 is about the cosmic mean at z = 14. The 
density fiuctuations may have effects on local Pop III 
formation by changing the LW background and metallic- 
ity level of the IGM. As discussed before, the strength of 
LW background during these redshifts might not change 
the Pop III formation significantly. In addition, because 
the mean distance between Pop III and halos is still much 
larger than the direct met al impact distance of a PISN of 
5 kpc (iWise et al.ll2012bl) . the metal enrichment by Pop 
III feedback crossing halos is still weak, except in halo 
clusters. In the case, the Pop III formation mostly de- 
pends on single halo properties, while the multiplicity of 
Pop III over halos is determined by the clustering of ha- 
los. Our results of Pop III distribution and multiplicity 
over halos should not change significantly in regions of 
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Figure 8. Complementary cumulative distributions of the metal- 
licity by volume (Top) and by mass (Bottom) inside the refined 
region at z = 15. The solid, dashed, dotted lines show the metal 
from Pop III stars, Pop II stars, and the total metal, respectively. 
The total metallicity distributions are still determined by the feed- 
back from Pop III stars at this epoch. 

different density fiuctuations. The population of Pop III 
is proportional to the halo density, so the number density 
from this simulation at z = 15 should be about 5 times 
that of the cosmic mean or close to the cosmic mean at 
z - 10. 

One important question about Pop III stars is when 
and in what environments their formation completely 
ends. Unfortunately, our simulation docs not show the 
cessation, or even the saturation of Pop III formation yet, 
though the increase of star formation rate slows down 
in the last tens of million years. The Pop III star and 
remnant number density is still low (lower than that of 
IWise et all (|2012bf ) 1 Mpc^ simulation) , and there still is 
a significant fraction of halos between 10^ and 10^ Mq 
that have no Pop III stars and are in a favorable LW 
radiation field. To show the sufficiency of metal poor gas 
inside the survey volume to continue form Pop III stars, 
we plot the volume and mass complementary cumula- 
tive distributions of the metallicity at z = 15 in Figure 
[Hi There are only about 0.65% of volume and 3.8% of 
the mass enriched above 10~^ Zq. For comparison, the 
IWise et al.l (l2012b[ ) 1 Mpc^ simulation has 0.32% of the 
volume and 2.6% of the mass at z = 10, and 0.89% of 
the volume and 7.6% of the mass at z = 7.28, respec- 
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tively, above 10~^ Z©. So it is reasonable to expect Pop 
III stars will continue to for m inside our s urvey v olume, 
but at a constant rate, since iWise et al] (l2012bO shows 
that the Pop III formation rate is constant from z = 12 
to 7. We are continuing the present simulation to lower 
redshifts to confirm or deny this expectation. 

This research was supported by National Science Foun- 
dation (NSF) grant AST-1109243 to MLN. JHW ac- 
knowledges support from NSF grant AST-1211626. The 
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covery Environment (XSEDE) by the National Institute 
for Computational Science, ORNL with XRAC alloca- 
tion MCA-TG98020N. 
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